Correcting ray optics at curved dielectric microresonator interfaces: 
Phase-space unification of Fresnel filtering and the Goos-Hanchen shift 
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We develop an amended ray optics description for reflection at the curved dielectric interfaces of optical mi- 
croresonators which improves the agreement with wave optics by about one order of magnitude. The corrections 
are separated into two contributions of similar magnitude, corresponding to ray displacement in independent 
quantum phase space directions, which can be identified with Fresnel filtering and the Goos-Hanchen shift, 
respectively. Hence we unify two effects which only have been studied separately in the past. 

PACS numbers: 05.45.Mt, 42.55.Sa 



C/2 

o 

• 1—1 

> 

Oh 

q 

O 
Oh 



(N 
> 

O 
O 

o 
o 

Oh 



Over the recent years it has become feasible to design opti- 
cal microresonators that confine photons by means of dielec- 
tric interfaces into a small spatial region not larger than a few 
micrometers 1 1, 2, 3]. Two promising lines of research are the 
amplification of photons by stimulated emission in active me- 
dia, which yields lasing action L L lAiSi. A iS i6. i7t i8i. i9t iT^ [Till . 
and the generation and trapping of single photons which can 
be used as carriers of quantum information [T?]. These ap- 
plications require integration of several components and in- 
terfacing with electronics, which are best realized in two- 
dimensional resonator geometries where the main in- and 
out-coupling directions are confined to a plane, and can 
be selected via the (asymmetric) resonator geometry. Fur- 
thermore, because of the requirements of mode selection, 
these applications favor microresonators of mesoscopic di- 
mensions, with size parameters kL ~ O(IOO) — 0(1000) 
(where L is the linear size, k — 2tt/\ is the wavenum- 
ber and A is the wavelength) which quickly puts these sys- 
tems out of the reach of numerical simulations. On the other 
hand, ray-optics predictions of the intricate resonator modes 
0, 13, '9', fll Eilil El Q can deviate substantially from 
experimental observations Q 0] and theoretical predictions 

aiiiniiiii. 

The purpose of this paper is to develop an amended ray 
optics (ARO) description which still idealizes beams as rays, 
but incorporates corrections of the origin and propagation di- 
rection of the reflected ray. We identify these corrections 
by utilizing quantum-phase space representations of the in- 
cident and reflected beam 11811 and relate them to the recently 
discovered Fresnel filtering effect fl^] and the long-known 
Goos-Hanchen shift |20]. These two effects have only been 
discussed separately in th e past (f or applications to microres- 
onators see, e.g., Refs. I5l ll l[l2l[l22il '). and their complemen- 
tary nature has not been realized. Moreover, their uniform 
analysis for all angles of incidence is known to pose consid- 
erable technical challenges In the phase- 
space representation, the Fresnel filtering and Goos-Hanchen 
corrections are simply determined by the position of maximal 
phase-space density. For the prototypical case of a Gaussian 
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FIG. 1: (Color online) Left panel: Gaussian beam reflected from a 
curved dielectric interface {kR = 100) separating regions of refrac- 
tive index tiq = 1 and n = 0.667. Light regions indicate high wave 
intensity. The angle of incidence xo = 42° is close to the critical 
angle x — 41.8°. Conventional ray optics predicts that the beam 
is specularly reflected at the point of incidence. In this paper we 
use phase-space representations to obtain a more accurate reflection 
law, which accounts for (i) the Goos-Hanchen shift (GHS) Acp of the 
reflection point along the interface and (ii) the increase Ax of the 
reflection angle due to Fresnel filtering. Both effects change the far- 
field radiation direction 9 (see the right panel, which exaggerates the 
corrections in order to clarify the notation). For the parameters in the 
left panel, ~ 7° and Ax ~ 1° (see Fig.|3}, resulting in a cor- 
rected ray which nicely reproduces the observed radiation direction. 



beam reflected from a constantly curved dielectric interface, 
we find that compared to conventional ray optics, the resulting 
ARO improves the agreement of the far-field radiation charac- 
teristics with wave optics by about one order of magnitude. 

The plan of this paper is as follows. We first quantify the 
corrections to conventional ray optics by phase-space repre- 
sentations of the wave in the near field of the interface. Then 
we incorporate these corrections into an ARO. Finally, we test 
the ARO by its predictions for the far-field radiation charac- 
teristics. 

Conventional ray optics assumes that beams have well- 
defined propagation directions and a precise point of impact 
on a sharp dielectric interface, and predicts that an incident 
beam is reflected specularly and locally at the interface i26ll . 
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In two dimensions, deviations from ray optics at curved inter- 
faces are apparent already at inspection of wave patterns such 
as shown in the left panel of Fig. [2 where the incident beam 
propagates from right to left. The wave pattern reveals that the 
reflected beam originates from a displaced position and prop- 
agates into a different direction than predicted by ray optics. 

We choose a coordinate system with origin at the center of 
the circle of curvature (see the right panel of Fig. [3. This cir- 
cle has the same radius of curvature R as the dielectric inter- 
face and touches it tangentially at the point of incidence. The 
incident beam propagates into negative x direction. For the 
comparison of wave optics to ray optics it is convenient to pa- 
rameterize the rays by Birkhoff coordinates (0, sin x), where 
(f) parameterizes the polar angle of the ray's intersection point 
with the interface, while R sin x is the impact parameter of the 
ray, where x is its angle of incidence. In this two-dimensional 
phase space, ray optics locates the incident and reflected rays 
in Fig.[2both at the same point (p = (po, sin x = sin xo, where 
furthermore (po = xo for the present case that the incident ray 
is oriented into negative x direction. 

In wave optics, the corresponding two-dimensional Gaus- 
sian beam is described by the wavefunction 

vl/i„(r,0)=^c„i/„(fcr)e""^ (1) 

rn 

_ g*(xo-0o-f )"i--^(m-fc-Rsin xo)^ (^-^ 

where are Hankel functions and w is the width of the 
beam in the polar angle (p. Since we are interested in the cor- 
rections in leading order of kR, we assume that the curvature 
is locally constant. Then the reflected beam has the wavefunc- 
tion 



(a) incident, near field 



(b) reflected, near field 



*rcfl(r,(/.)-5]c+i/+(fcr)e™^ 



(3) 



^+ ^ ^- nj:^{nkR)H;^{kR) - H'^ {kR)J^{nkK) 
^" ^" H+'{kR)Jm{nkR) - nJ^{nkR)H+{kR) ' 

where J™ denotes the Bessel function and n is the refractive 
index on the other side of the interface. 

In order to analyze the exact wave pattern in the phase space 
of Birkhoff coordinates we utilize Husimi representations in 
the near field of the dielectric interface. These Husimi func- 
tions are obtained by overlapping the incoming and reflected 
beam at the interface with a minimum uncertainty wave packet 
centered around {(p, sinx) Ll&l . 



, sinx) = cosx X 



i{m—kR sin x)4^~ [m — kR sin x) 



(5) 



The Husimi phase-space representations of the wave pat- 
tern of Fig.[2is shown in Fig.|2 The width w = yJir/kR of 
the incident Gaussian beam is chosen such that it yields an op- 
timal approximation of a classical ray with comparable uncer- 
tainties in the propagation direction and the point of impact. 




FIG. 2: (Color online) Near-field Husimi phase space representa- 
tions js} of the wave pattern of Fig. Q The left panel shows the 
Husimi function 7i~ of the incident beam, while the right panel 
shows the Husimi function Ti"*" of the reflected beam. The crosses 
+ indicate the ray-optics prediction for the point of highest phase- 
space density, which is accurate for the incident beam, but not for 
the reflected beam. The displacement into (p direction can be related 
the Goos-Hanchen shift, while the displacement into sin x direction 
is the consequence of Fresnel filtering. 



This results in the almost-circular phase-space density in the 
left panel. The location of the maximal phase-space density 
corresponds well with the ray-optics prediction (xcsinxo)^ 
indicated by the cross +. The phase-space representation of 
the reflected beam is shown in the right panel. Clearly the po- 
sition ((/)rnax, sin Xmax) of the maximal phase-space weight of 
the reflected beam is displaced from the ray-optics prediction, 
as had to be expected from the inspection of the wave pattern 
in Fig.n 

The displacement in sinx direction can be explained by 
Fresnel filtering, which was introduced by Tureci and Stone 
lIsl fToil : A realistic beam has an uncertainty > 1/ (kRw) of its 
propagation direction which results in a spreading of the angle 
of incidence. The angle of incidence is further spread because 
of the curvature of the interface over the focal width Rw. The 
Fresnel reflection coefficient displays an angular dependence 
which favors the reflection of wave components with a larger 
angle of incidence. This increases the beam's angle of reflec- 
tion, by an amount which we identify with the displacement 



Ax = Xn 



Xo- 



(6) 



The displacement into (p direction can be interpreted as a 
Goos-Hanchen shift (GHS), first discovered for planar inter- 
faces in 1947 ^ (for recent works see Refs. IlllE^Eilli 
I30I1). This shift arises from the penetration of the evanescent 
wave into the optically thinner medium 1.23, 24„1 . We identify 
the resulting lateral displacement of the reflection point along 
the physical interface with 



A(P ■ 



(7) 



The angle-of-incidence dependence of A(p and Ax is 
shown in Fig.|3] Both corrections are most pronounced around 
the critical angle of incidence x' ~ 41.8°, and are sizeable 
effects even for rather large values of kR. Beyond the criti- 
cal angle, approaches the classical result for the GHS by 
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(a) incident, far field 



(b) reflected, far field 



kR= 100 



kR=400 



30° 



40° 



50° 



60° 



FIG. 3: (Color online) Angle-of-incidence dependence of the Goos- 
Hanchen shift and the Fresnel filtering correction Ax in the near 
field of a curved interface with kR =100 (top panel) and kR — 400 
(bottom panel). The remaining parameters are as in Fig. Q The 
black line shows the classical result for the Goos-Hanchen shift by 
Artmann Isill. 



Artmann 1123113 111 , which is derived in the regime of total re- 
flection X > x' at a planar interface. 

ARO consists in propagation of the reflected beam with 
point of reflection given by 0^^° — _|_ ^jjjj angle of re- 
flection given by x^^*-* — Xo + ^Xj resulting in a propagation 
direction 



qARO 



= (/-o + Xo + A(/. + Ax 



(8) 



(see the right panel of Fig. [Q. Note that the corrections 
and Ax both have been determined in the near field of the 
interface [see Eqs. (I6I7H . Hence, within the idealization of 
beams by rays, ARO agrees exactly with wave optics in the 
near field of the interface. The question is then whether the 
ARO ray parameters deliver accurate predictions also in the 
far field (where the beam may encounter another optical com- 
ponent or a detector). Hence, we now test the accuracy of 
ARO by examination of its predictions for the far-field radia- 
tion characteristics. 

Figure 0] assesses these predictions for the wave pattern of 
Fig.[2by means of Husimi phase-space representations of the 
incident and reflected beam in the far field of the interface. 



,sinx) = 



E 



± i{7n~kB. sin x)0T''''^^" /2— {m — kR sin x)'^ 



(9) 



In the far field, the phase-space coordinate (j) coincides with 
the propagation direction 9, while sinx is still related to 




FIG. 4: (Color online) Same as Fig.|2| but for the far field, where 
the Husimi representations are given by Eq. J9}. The diagonal 
cross X indicates the ARO prediction for the reflected beam, which 
incorporates the Goos-Hanchen shift and Fresnel filtering. 



the impact parameter i?sinx (this coordinate is preserved 
because of angular-momentum conservation with respect to 
the center of the circle of curvature). The incident beam 
propagating into negative-x direction is thus represented by 
phase space coordinates (</), sinx) = (O,sinxo)- Ray-optics 
predicts that the reflected beam has phase-space coordinates 
((/),sinx) — ,smxo), where 6*^*-* — <j>o + Xo, while 
ARO predicts that the reflected beam is located at {(f), sin x) = 
^^ARO gjj^^ARO-)^ The position (^max, sinxmax) of the max- 
imal phase-space density of the reflected beam in the far field 
(right panel of Fig. ^ indeed corresponds well to the ARO 
prediction (indicated by x), but deviates distinctively from 
the ray-optics prediction (indicated by +). 

In Fig. |5] the far-field radiation direction 9 is analyzed as a 
function of the angle of incidence. One of the curves is the 
deviation A9^'^ = i^max — d^'^ of the observed radiation 
direction from the prediction of conventional ray optics. For 
kR = 100, the maximal deviation is « 12.5° and occurs about 
4° below the critical angle of incidence. The plot also shows 
the deviation of ARO, A9^^° = A9^^ -A(j)-Ax ■ Itis seen 
that ARO improves the agreement to 2° close to the critical 
angle and agrees even better away from it. For larger size 
parameters kR — 400, the maximal disagreement between 
ray optics and wave optics drops to « 3.5° and occurs at about 
1° below the critical angle of incidence. The ARO prediction 
agrees within 0.3° around the critical angle, and the agreement 
is almost perfect away from it. 

To summarize, we developed a systematically amended 
ray-optics description of the reflection of Gaussian beams 
from the curved dielectric interfaces of optical microres- 
onators. This description incorporates the Goos-Hanchen shift 
of the reflection point along the interface and the Fresnel- 
filtering enhancement of the angle of reflection. The cor- 
rections were determined by analysis of exact wave-optical 
beams in a phase space where one coordinate is associated 
with the point of incidence or reflection along the interface, 
while the other one is related to the angle of incidence or re- 
flection, respectively. Fresnel filtering and the Goos-Hanchen 
effect displace the reflected beam along independent phase- 
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FIG. 5: (Color online) Angle-of-incidence dependence of the devia- 
tion of the far-field radiation direction 9 from the predictions of ray 
optics (RO) and amended ray optics (ARO). Top panel: kR = 100. 
Bottom panel: kR = 400. The remaining parameters are as specified 
in Fig.Q] 



space directions. Hence, these displacements in principle ex- 
haust all possibilities of amending ray optics while still keep- 
ing the basic assumption of propagation along straight lines in 
optically homogeneous media. 

Amended ray optics is applicable to microresonators with 
smooth boundaries where the dimensionless radius of cur- 
vature kR is large, which is realized in most experiments. 
This includes th e popul ar examples of multipole deforma- 
tions |4, 5, Hill Hill El or stadium geometries 
Complementary techniques exist to describe the diffraction of 
beams at sharp corners where formally kR = Is^l : see Ref. 

for an application to hexagonally shaped resonators 0]. 
It remains to be seen whether both techniques can be inter- 
laced to describe geometries which combine both curved in- 
terfaces and sharp corners lliol [Till ; moreover, whether both 
techniques can be unified in the challenging regime of a local 
curvature with kR = 0(1). 
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